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An Infrared Nebula Associated with 6 Cephei: 
Evidence of Mass Loss? 
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ABSTRACT 

We present the discovery of an infrared nebula around the Cepheid prototype 
5 Cephei and its hot companion HD 213307. Large scale (~ 2.1 • 10^ AU) nebu- 
losity is detected at 5.8, 8.0, 24 and 70 /im. Surrounding the two stars, the 5.8 
and 8.0 /im emission is largely attributable to Polycyclic Aromatic Hydrocarbon 
(PAH) emission swept from the ISM by a wind originating from 6 Cephei and/or 
its companion. Stochastically heated small dust grains are the most likely source 
of the 24 and 70 fim extended emission. The 70 /xm emission, in particular, re- 
sembles a bow shock aligned in the direction of the proper motion of S Cephei. 
This discovery supports the hypothesis that 6 Cephei may be currently losing 
mass, at a rate in the range ~ 5 • 10^^ to 6 • 10^^ M© yr~^. 
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Introduction 



5 Cephei is a remarkable s tar. After John Goodricke recognized its variability in 1784 
( IGoodricke &: Staatsoperlll786l ). 5 Cephei has lent its name to one of the most important 
classes of variable s tars: Cepheid s. Since the discovery of the period-luminosity (PL) relation 
(or "Leavitt law", iLeavittI llQOSi ). Cepheids and their prototype 5 Cephei have assumed a 
fundamental role as the "first rung" of the extragalactic distance scale. Cepheids are also 
a crucial benchmark for evolutionary models of intermediate-mass stars in the He-burning 
phase. Despite their importance, however, there are still a number of outstanding puzzles in 
the theoretical understanding of the stars pulsating like 5 Cephei. 

The most important iss ue is the d iscrep a ncy between t he theoretical and dynamical mas s 
of Cepheids. First noted by Christyl ( 1968 ). Stobie ( 1969 ) and later by Fricke et al. ( 1972 ) 



this issue has been pa rtially solved 



of radiative opacities ( ISeaton et al 



IMos ^alik et al. 19921) with th e adoption of modern sets 



19941 [Rogers fc IglesiasI Il992l ) . However, several recent 



investigations : 


ocused 


Beaulieu et al. 


2001 





Keller fc Wood 2006; Evans et al. 



Bono et al 



2008h 



2001, 2002 



suggest 



that s uch a discrepancy still amounts to 10-15%. For example, according to ICaputo et al. 
( I2OO5I ) the pulsational mass of 5 Cephei is Mp ~ 4.5 M0, while the evolutionary mass of 



the star, based on the period-color-luminosity relations derived by the same authors, can 
be as high as Mev — 5.7 Mq. One possible solution for this problem is mas s loss over the 



lifetime of these stars. Even though com monly used semi-empirical relations (IReimers et al. 



I975I : iDeJager &: NieuwenhuijzenI 119971 ) do not predict enough mass loss to account for the 
Cepheids' missing ma ss, these same relations are clearly inadequate to describe the mass loss 



( lYong et al.ll2000l : ICastellani et al 



i n red giants (see e.g. IWillsonI 12000) and in globular cluster stars on the Horizontal Branch 



20061 ). These relations may similarly underestimate the 



mass loss rates in intermediate-mass stars like Cepheids. 

Searches for evidence of current or past mass loss by Cepheids and their progenitors 
have been co nducted at different ends of the electromagnetic spectrum. In the ultraviolet, 
Deasyl (119881 ) observed a large sample of Cepheids with the International Ultraviolet Explorer 
(lUE) satellite, fiiiding m ass loss rates ranging from 10^^'^ to 10"*" M(n yr~^. In the far- IR, 
McAlary &: Welch! (119861 ) found evidence in the IRAS photometry (IBeichman et al.lll988[ ) of 
very cool dust (T^ < 50 K) around two cl assical Cepheids (RS Pup and SU Cas) associated 
with reflection nebulosity. More recently, iNeilson et al.l (|2009[ ) measured mass loss rates of 
10~^^-10~^ Mq yr~^ in a sample of Large Magellanic Cloud Cepheids in the mid-IR. In the 
near-IR, K band i nterferometric obseryations detected circumstell a r she l ls aro und nearby 



classical Cepheids (iMerand et all 120061 l2007t iKervella et aP bood . l2008l . l2009h . including 



5 Cephei. A difference in the absolute distance of 5 Cephei when moving from V, I to V, K 
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bands was also found by (INatale et al.ll2008l ) by performing a detailed fit between predicted 
and observed light, radius, and radial velocity curves. 

These detections offer tantalizing evidence in favor of Cepheid mass loss, but are not 
conclusive. In the case of the far-IR excess, the large beam of IRAS (~ 5') does not allow 
separation of local dust emission from background "Galactic cirrus" . The near-IR shell s 
detected with interferometers have radii as small as ~ 2 stellar radii (iMerand et al.ll2006l ). 
and their origin is not well understood. To resolve the se issues, we have obtained Spitzer 
Space Telescope ( Werner et al. 2004a : Gehrz et al. 20071) observations o f 29 Galactic Classi- 
cal Cepheids with the Infrared Array Ca niera (IRAC, iFazio et al.ll2004l ) and the Multiband 
Infrared Photometer for Spitzer (M I PS. iRieke et al.l 120041 ). The results of this investiga- 
tion are presented i n iMarengo et all J201of r (PL and search for mid-IR color excess) and in 



Barmby et al.ll201ll (search for extended emission). 



In this paper we want to focus on 6 Cephei. With a distance of 0.273±0.011 kpc 
(IBenedict et al.l |2002| ) and a fundamental mode period of 5.37 d, 5 Cephei is the second 



nearest Cepheid to the Solar System (only Polaris, a first overtone Cepheid, is closer). 
S Cephei is also known to be part of a multiple star system. The B7-B8 II I- V main sequence 



star HP 213307, long suspecte d to be a companion of 6 Cephei (IFernie 



Vitrichenko fc Tsarevskii 



1982 



Benedict et al.ll2002l). Located 4 0^^ 



19691 ) shares proper motion with the Cepheid (IHoffleit &: Jaschek 



966 



Worle 



1966 



11, 000 AU) south of 6 Cephei, HD 213307 may 
itself be a binary star (iBenedict et al.ll2002l ). with a low mass companion (~ 1.6 Mq, FO V) 



on a circular orbit with a radius of approximately 7 ra as (1.9 AU). The three stars are part 
of the moving group Cep 0B6 (Ide Zeeuw et al.lll999l ). together with another 18 members 
that are possibly part of an old OB association. The earliest spectral type (B5 III) among 
the stars in the association suggests an age of ~ 50 Myr, consistent with the es timated age 
of S C ephei {tage — 60 Myr according to the period-age- color relation derived by iBono et al. 



20051 . and using 6 Cephei colors from lStorm et al.ll2004l ). The proximity of S Cephei to the 



Sun, and the presence of a hot companion which can serve as a "probe" of the circumstel- 
lar environment of the system, make 6 Cephei the ideal test case to search for evidence of 
Cepheids mass loss. 

We report on the discovery of extended circumstellar emission around 5 Cephei and 
its companions, detected at IRAC and MIPS wavelengths. This emission appears to be 
physically associated with 6 Cephei, and may be evidence of active mass loss from this star. 
We detect a 70 /im arched structure with its axis aligned with the relative proper motion of 
6 Cephei with respect to the interstellar medium (ISM), and diffuse and filamentary emission 
at 5.8, 8.0 and 24 /im. In section [2] we present our Spitzer observations, while in section [3l 
m and [5] we discuss a possible explanation for the origin and the morphology of the observed 
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structures, and estimate the mass loss rate required for their formation. In section [6] we 
summarize our resuhs and conclusions. 



2. Observations and Data Reduction 



S Cephei was observed with IRAC and MIPS in Cycle-3, as part of our General Ob- 
server program with PID 30666, and in Cycle-5 (IRAC only) as part of our Guarantee 
Time Observer program PID 50350. The GTO observations were meant to provide a second 
photometric epoch, and deeper observations aimed at being sensitive to faint circumstellar 
emission. The observations were executed on 2006 August 10 (JD = 2453957.601) and 2008 
September 23 (JD = 2454732.832), using the IRAC full frame mode. The 2006 observations 
were obtained with a 5 point Gaussian dither pattern in IRAC "stellar mode", with total 
integration times of 1.0 sec (3.6 and 4.5 /xm) and 6.0 sec (5.8 and 8.0 fim). The 2008 data 
were instead acquired with a 36 point Reuleaux dither pattern with individual exposures 
of 1.6 sec integration time each, for a total exposure of 43.2 sec in each band. The MIPS 
observations were obtained on 2007 July 19 (JD = 2453935.518) in the 24 and 70 fim bands 
with the Photometry Astronomical Observation Template. The total on-source time was 
48.2 sec and 37.7 sec at 24 fim and 70 /im respectively. 



2.1. IRAC images and photometry 



The IRAC data were reduced starting from the Basic Calibrated Data (BCDs) generated 
by the Spitzer pipeline versions S14.4.0 and S18.7.0 (2006 and 2008 datasets respectively). 
Mosaic in aages with a pixel sca le of 0.8627"/pix were created using the IRACproc post-BCD 
software (ISchuster et al.ll2006[ ). 



In all IRAC images 6 Cephei is sa turated and we have d etermined its photometry using 
the PSF- fitting technique described in lMarengo et al.l ( l2010l ). The same technique was used 
to estimate the photometry of HD 213307 in the 2008 dataset, because of column pulldown 
(3.6 and 4.5 /im) and banding (5.8 and 8.0 fim) from the primary star overlapping with the 
position of the companion, making it impossible to derive accurate aperture photometry. For 
the 2006 data the HD 213307 brightness was determined with standard aperture photometry. 
The IRAC magnitudes so derived are listed in Table [1] 

The IRAC images are shown in Figuredl with 6 Cephei and HD 213307 PSF-subtracted. 
For the first time, these stars are found to be surrounded by extended infrared emission, 
brightest at 8.0 /im and marginally detected at 5.8 /im. The emission is diffuse, mostly 
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contained within the dashed arc plotted in the figure, and brighter in the region between the 
two stars. The extended infrared emission is on much larger spatial scale (up to ~ 1.3', or 
~ 2.1 ■ 10'^ AU at the 5 Cephei distan ce) than the ci r cums tellar near-infrared shell detected 



in the interferometric observations by lMerand et al.l (120061 ) with a ~ 2 R^, ~ 0.4 AU radius 



The surface brightness of the extended emission, measured in the deeper second epoch 
in the four boxes shown in Figure [H is listed in Table [2l The sky background level has 
been measured in the two dashed boxes also shown in Figure [H To take into account the 
background offset induced by array muxstriping, we have used sky box si for apertures 1 
and 2, and sky box s2 for apertures 3 and 4. We have corrected the IRAC surface brightness 
with their extended source aperture corrections of 0.91, 0.94, 0.70 and 0.74 at 3.6, 4.5, 5.8 
and 8.0 fim respectively. Where no extended emission was detected, we listed the 3a limits. 



2.2. MIPS images and photometry 

We reduced th e MIPS data using the MIPS instrument team Data Analysis Tool 



(iGordon et al.ll2005l ). The processing of the 24 /xm data was straightforward, and the fi- 
nal image is shown in Figure [H To remove column-like instrument artifacts, the 70 fim data 
needs to be time filtered with the source region excluded: this operation is complicated in 
our images because of the narrow field of view, and extended luminosity covering a substan- 
tial fraction of the detector. In Figure [2] we show the resulting filtered images using three 
different filtering choices. Panel a shows the mosaic without any filtering applied (leaving 
the strip pattern noise). The center and right panels show instead the result of excluding 
the default source point region size (6, similar to the default region chosen by the Spitzer 
Science Center post-BCD pipeline) and excluding the whole nebulosity area (panel c). 

The 24 /im nebulosity is morphologically similar to the one in the IRAC 8.0 /im map, 
and is also mostly contained within the arched envelope plotted in Figure [H even though 
there are faint extended lanes at larger spatial scales. It is also brightest close to 6 Cephei 
and its companion, showing the same structures detected at 8.0 /xm. 

The 70 /im image shows strong extended emission on the N-E side of S Cephei, stronger 
along the arched envelope marked in Figure [H The detailed morphology of this structure 
depends on the choice for the filtering exclusion region. The default choice highlights the 
high spatial frequencies allowing us to follow the arched structure through the whole width of 
the array. The larger exclusion region instead enhances the low spatial frequencies, showing 
better the faint nebulosity in between the star and the cusps of the arc. The actual shape of 
the nebulosity is most likely a combination of the two, but a 70 /xm image with larger field of 
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view would be required to properly filter the map and assess its background emission level. 
The 24 /xm and 70 /im magnitud es of 6 Cephei and its c ompanion, derived with standard 



aperture photometry as described in iMarengo et al.l (120101 ). are listed in Table [H At 70 fim 



6 Cephei is detected with less than 3a significance (93 ± 50.5 mJy, from the default filtered 
image, which best preserves the point source). Within the uncertainty, the color of the star 
(with respect to the simultaneous 24 /xm band) is zero, indicating absenc e of significant 



amou nt of warm and cold dust near the star (confirming the results given in iMarengo et al. 



20101). 



We have estimated the 24 and 70 fim surface brightness of the extended emission in 
the same four boxes shown in Figure [H The background level was estimated by averaging 
the surface brightness values in the two sky boxes, also marked in Figure [H The surface 
brightness is listed in Table [H For the 70 /xm emission, we list the two values corresponding 
to the different filtering options described above. 



3. The 6 Cephei Nebula 



With a Galactic latitude of +0.53°, 6 Cephei is directly projected against the outer 
Galactic plane, in an area with diffuse infrared-emitting background "cirrus". Inspection 
of wide field 24 and 70 /im part of the ''Spitzer Mapping of the Outer Galaxy" survey 
(SMOG), from the Spitzer archive, shows that the diffuse emission centered around S Cephei 
is contiguous to a system of larger Galactic structures. Understanding if the nebulosity 
around 6 Cephei is physically associated with the star, rather than a chance superposition 
with background interstellar clouds, is critical to investigate the mass loss history of the 
star. As part of the Spitzer Cycle-5 p rogram PID 50346, we proposed to obtain a Infrared 
Spectrograph (IRS, iHouck et al.ll2004l ) spectral energy distribution of the 6 Cephei extended 
emission in key areas around the star, to investigate the presence of diagnostic features within 
the IRAC and MIPS bands, and thus identify the nature of the circumstellar matter. The 
exhaustion of the spacecraft LHe refrigerant at the end of the cryogenic mission, however, 
prevented this part of the program from being executed. Program PID 40968 on 2007 
December 10 obtained an IRS high resolution spectrum centered on 6 Cephei. The data 
however do not show other emission than the strong photospheric flux. 

In absence of spectra for the extended emission, we compared our 6 Cephei IRAC and 
MIPS ima ges with the observat ions of a similar source, for which IRS spectra are available. 
Following iKervella et al.l (120091 ) . we adopted the dusty reflection nebula NGC 7023 as an 
"analog" for a Cepheid precursor surrounded by circumstellar nebulosity. The central source 
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in this system, the binary star HD 200775, has a B2Ve composite spectral type, similar 



to th e 6 Cephei B7-8 companion HD 213307. With masses of 11 and 9 Mq ( lAlecian et al. 



20081 ). the two stars at the center of NGC 7023 may become Cepheids themselves within a 
few tens of million years. Studying the diffuse emission in NGC 7023 thus offers the chance 
to probe a system both representative of an early version of 6 Cephei, and powered by a 
central stellar engine similar to the S Cephei companion. 

NGC 7023 was observed with all Spitzer instruments as part of the Early Release Ob- 
servation programs 717 and 1093. The central stars are embedded in a reflection nebula 
characterized by bright filaments surrounding an hourglass-shaped cavity (filled with CO 
emission) probably formed by outflow activity in the HD 200775 past. The IRS spectral 
energy distr ibution in a location along one of the bright filaments has been presented and 
discussed in IWerner et al.l (l2004b( ). We have derived the surface brightness of the filaments 
in all IRAC and MI PS images in the same filament location (position "A" in Figure 2 of 
Werner et al.ll2004bl ). and listed them in Table |2l This location is situated at a distance of 
~ 2.0 ■ 10^ AU from th e central stars (using a distance of 430 pc for NGC 7023, according to 
van den Anckerl 119971 ). similar to the radius of the bright 70 fim arc in the S Cephei system. 
The availability of the complete infrared spectrum at this location allows us to determine 
the prevalent mechanisms of emission responsible for the NGC 7023 diffuse luminosity in the 
IRAC and MIPS bands and, by extension, provide diagnostic tools for the unknown emission 
around 6 Cephei. 



Werner et al.l (l2004bl ) showed that the bulk of the flux in the IRAC bands at 5.8 and 



8.0 /zm is due to strong PAH and other aromatic emission features. Fainter emission detected 
at 3.6 and 4.5 /im is likely related to weaker H2 lines. The stronger flux in the MIPS bands, 
in contrast, is associated with thermal emission from Very Small dust Grains (VSGs) that 
are stochastically heated by single UV photons generated by the central B stars. This 
interpretation is confirmed by comparing the S'8.0/5'24 and 82^ / Sjo surface brightn ess ratio 
of the filament with models of interstellar dust computed by iDraine fc Lil (120071 ). These 
models cannot be strictly applied because of the different assumption in the spectrum of the 
input energy field: two B stars in the case of the N GC 7023 nebu l a and a scale d interstellar 
radiat ion field derived for the solar neighborhood by lMathis et al.l (119831 ) for the lDraine fc Li 
( l2007l ) models. It is nevertheless instructive to note that the flux ratios in the NGC 7023 
filament can be reproduced by assuming an irradiation of the nebular material by a field 
equal to ~ 100 times the interstellar radiation field in the solar neighborhood, and a PAH 
dust fraction g ~ 3%. 



The surface brightness ratios of the NGC 7023 filament, and of the four locations in the 
5 Cephei nebula listed in Table |2] are plotted in Figure [31 The figure shows that the 5 Cephei 
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emission has a Sf^ o/S^i significantly smaller than NGC 7023, suggesting a PAH content 
q < 1% (based on lDraine fc Lil 120071 ). The 6 Cephei nebula S2aISiq ratio along the 70 /im 
arc is as high as ~ 3 times the corresponding ratio for the NGC 7023 filament, suggesting a 
temperature distribution of the dust grains around the Cepheid peaked at higher values. It 
is worth noting that the S24,/S7o flux ratio (and hence the dust temperature) appears to be 
loosely correlated with the distance of each box from the two stars in the S Cephei system. In 
particular, the flux ratio in box 1, located on the bright filament in between the two stars, is 
more than 4 times larger than in the other three boxes: this is strong evidence supporting the 
hypothesis that the main energy input for the 6 Cephei nebula is the two central stars, and 
thus that the nebula itself is local, and not a background cloud in a chance alignment. We 
do not detect any emission at 3.6 and 4.5 fim. Nevertheless, we cannot exclude the presence 
of shocked H2 molecules around 6 Cephei, since H2 line emission would be undetected in 
our images if its S^ q/Ss q and S^^^/Ss.o ratios were similar to what we have measured in the 
NGC 7023 filament. 



A far-IR Bow Shock? 



The envelope of the 70 /im emission is shaped like a paraboloid centered on 5 Cephei 
and its companion. Structures of this kind are often associated with bow sh ocks that form 



due to the interaction of a stellar wind with the surrounding ISM (see e.g. IWareing et al. 



20061 ). Bow shocks can be bright in the thermal infrared, when the accumulated interstellar 
material, or the stellar wind itself, is rich in dust. Examples of these structures around 
evolved stars abound, including the cas e of the Asymptotic Giant B ranch (AGB) star R Hya 
(jUeta et all I2OO6I : IWareing et all 120061 ) and the supergiant a Ori (jUeta et al.l |2008aJ) . The 
Mira variable class prototype o Ceti also shows a complex infrared bow shock structure, 
associated with a UV and infrared-bright "coraetary t ail" that is generated as the star plows 
at high velocity through its local ISM ( lUetal l2008b| ). If the far-IR structure surrounding 
S Cephei is indeed a bow shock, its axis is expected to be aligned with the direction of the 
star's motion through the ISM. 



Following the prescription of IJohnson fc SoderblomI ( 1l987l ). we have computed the com 



ponents of the Galactic peculiar space motion of the two stars, {U, V, W)pec- For these 
calculations we assumed the following parameters f or 6 Ce p hei (H P 213307): heliocentric 
radial velocity —16.8 km s~^ (—21.0 km s~^) from IWilsonI ( 119531 ): proper motion in right 
ascension 16.47 mas yr^^ (16.48 mas yr^^) and proper motion in decl ination 3.55 mas yr~^ 
(4.70 mas yr^^) from the Hippar cos database (IPerrvnaan et al.lll997| ). We used the paral- 
lax of 3.66 mas (3.65 mas) from Benedict et al.l ( 2002 ). With the exception of the radial 
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velocity, the values for the two stars are al l consistent to within unc ertainties. Correction 
for the solar motion using the constants of iDehnen &: Binneyl (Il998l ) yields {U,V,W)pec = 
(-6, -16, 0) km s-i for 5 Cephei and {U, V, W)pec = (-5, -21, 1) km s"^ for HD 213307. Pro- 
jecting back into an equatorial reference frame then gives {Vr,a,6)pcc = ( — 14,8,5) km s~^ 
and (Vr, a, 6)pec = ( — 18, 8, 6) km s~\ respectively. This implies a space velocity for 6 Cephei 
of V^pace ~ 17 km s~^ along a position angle of 59° and a space velocity for HD 213307 of 

space ~ 



21 km s ^ along a position angle of 52°. 



Figure [T] shows that the velocity vector of 6 Cephei with respect to the motion of its 
local ISM, is roughly aligned with the symmetry axis of the 70 fim paraboloid emission. 
While this is not a proof that the shape is indeed the result of a bow shock, it strongly 
suggests that this may be the case. If a bow shock forms in front of 6 Cephei, conservation 
of momentum implies that the star must be losing mass, in the form of a wind interacting 
with the ISM. 

How much dust mass is associated with the 70 /im emission? If we assume that the 
dust is optically th in, we can es t imate the total mass of the dust responsible for the 70 fim 
emission, following lEvans et al.l ( l2003l ): 



4.97 X 10-^Vz. 



KyBy{Td 



where is the total flux of the extended emission at 70 /xm in J y, D is the distance in kpc, 
= 56 cm^ g~^ is the dust opacity from lOssenkopf fc Henning] (119941 ) for a Mathis-Rumpl- 



Nordsieck dust distribution (IMathis et al 



19771 ). By{Td) is the Planck function in c.g.s units 



and Td is the temperature of the dust. 

Even though, as noted in Section [3l the small dust grains responsible for the far-IR 
emission are most likely stochastically heated, we can still estimate the mode of the grain 
temperature distribution as the color temperature of the S'24/ 5*70 flux ratio. Boxes 2, 3 and 4, 
lying on the 70 /im bright arc, have color temperature in the range of 85-110 K, comparable 
with the color temperature of the NGC 7023 filament (~ 80 K). The color temperature of 
the dust in box 1, placed in between 5 Cephei and HD 213307, is however much higher 
(~ 180 K), as expected if heated by the radiation from the two close- by stars. 

The total 70 //m flux density of the arched emission, measured in the dotted rectangular 
annulus in Figured] is between 1.0 and 2.8 Jy, depending on which data reduction product 
is used. We used a rectangular aperture to exclude the flux from the marginally detected 
star. From Equation [H using the average color temperature Td ~ 100 K along the arched 
structure, and fu ^2 Jy, we have Md ~ 6 • 10"'' Mq. If we assume a gas to dust mass ratio 
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of ~ 100, typical of ISM and circumstellar dust, we obtain a total mass of gas and dust of 
^ 6 • 10~^ Mq. a minimum age of this structure is given by the time required by a stellar 
wind to reach the arc from 6 Cephei (or its companion). In the case of 6 C ephei, the wind 



speed is expected to be comparable to the escape velocity (~ 100 km s ^; e.g. IWelch fc Duric 
19881 ). i.e. several times higher than is typical for AGB stars. This implies a minimum age 



of > 10^ yr for the structure, and a maximum mass loss rate of M ^ 6- 10 ^ M , ^ yr ^. These 



value is at the upper end of the Cepheids mass loss rates proposed by iDeasyl (119881 ). 



We can also derive a minimum mass loss rate required to produce the obse rved bow shock 
struc ture using the standard stellar wind/ram pressure balance condition (e.g.. lRaga fc Canto 
2008h : 



M 



(2) 



Here dg is the standoff distance of the bow shock, po is the ambient ISM density, is 
the space velocity of the star relative to the ambient medium, and is the wind outflow 
velocity. From our Spitzer images we measure dg ^ 1.3 arcmi n (0.103 pc). We assume t hat 
the surrounding medium is pure atomic hydrogen and, from iDickey &: LockmanI (jl990l ). a 
typical particle density in the Galactic plane is ~ 0.55 cm~'^, implying po ~ 9.2-10"^^ g cm~^. 
Finally, adopting an expected outflow velocity for the Cepheid wind of ~ 100 km s~^ 
(see above), we find M ^ 5 ■ 10~^ Mq yr~^. While this value can be considered only 
approximate given the uncertainties in the wind sp eed for the st ar and the local ISM density, 
it is nonetheless consistent with the predictions of iDeasyi (1l988l ) for Cephe id mass-loss rates 
(albe it significantly higher than the Mpuu — 3.6-10"^" M0 yr~^ predicted by lNeilson fc Lester 
20081 for this star). In particular, if 6 Cephei sits in a over-dense area of the Galactic plane, 
as the SMOG and IRAS maps imply, the ISM density (and thus the minimum M) could be 
higher. 

The fact that the mass-loss rate derived using this approach is one magnitude lower than 
the estimate computed above, using the observed properties of the dust emission, suggests 
that part of the bow shock may be comprised of swept up interstellar material rather than 
entirely shed from the star. Our images in fact do not allow us to determine how much of the 
detected mass is coming from the star, rather than being swept from the ISM by the stellar 
wind (see also discussion below on why at least the PAH emission detected at 8.0 pm is 
likely of interstellar origin). The presence of a bow-shock structure, however, argues in favor 
of a stellar wind pushing against the ISM, and thus at least part of this mass ought to come 
from the star, with a mass loss rate of at least 5 ■ 10^^ Mq yr~^ and up to 6 • 10^*^ Mq yr^^. 
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Given the high luminosity of the Cepheid and its companion, it is fair to ask if the 
observed dust-free bubble around the stars may be generated by effect of radiation pressure 
alone, without the need of invoking a Cepheid winds. The radiation pressure has an inverse 
square law variation and therefore acts like a repulsive "Coulomb" force. If an interstellar 
grain of mass m has speed v = 17 km s~^ relative to the star, then the distance of closest 
approach D is given by: 



—mv^ = — - — Tca^Q (3) 

where is the luminosity of the Cepheid (~ 2, 000 Lq), c is the speed of light, a is the radius 
of a grain and Q is a dimensionless number to express the effective momentum coupling to 
the geometric cross section. For a spherical grain of radius a and density — 3 g cm~^ 
(typical value of astronomical silicates): 



D = (4) 

Given that dust grains are opaque at optical and near-IR radiation, we can assume 
Q = 1. According to equation H] a distance of closest approach ~ 0.1 pc is obtained for 
dust grains with radius a ~ 0.11 fim. 

However, based on the discussion in Section [31 the 24 and 70 fim emission is most likely 
due to stochastically heated VSGs, that have a typical radius of ~ 10~^ /xm. Such grains 
would be blown out by radiation pressure at a distance two order of magnitudes larger than 
the observed D ~ 0.103 pc. Furthermore, the opening angle of a dust-free bubble created by 
radiation pressure would be much larger than the one observed: for the direction orthogonal 
to the space velocity of the star with respect to the ISM, the dust grains will have a zero 
radial component velocity, and will be blown out at a distance much larger than D derived 
from Equation HI in contrast to the narrow arc observed in Figure [1] The small size of the 
emitting grains, and the shape of the arched structure, imply that gas drag is likely to play 
a strong role in preventing the dust grains from being blown out by the radiation pressure. 
While radiation pressure is likely to play some role in the dynamic of the system, we thus 
still favor the hypothesis that the arched structure is due to a mechanical force acting on 
both interstellar gas and dust, due to a wind originating from the stars in the system. 

Detection of emission lines from the shocked material along the arc would confirm that 
the structure is indeed a bow shock generated by a stellar wind pushing on the ISM. The 
power P in the stellar wind is given by |Mf^. With ~ 100 km s~^ and a mass loss rate 
in the range of 5 • 10"^ to 6 ■ 10"*^ M© yr^^ P 4 ■ 10"^ to 5 • 10"^ jhis energy will 
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likely be dissipated in a radiative shock, and emission lines (among which H2 or CO) could 
be observable. This was one of the main motivations for our Spitzer /1R.S cycle-5 program, 
that was unfortunately not executed before the end of the cryogenic mission. We have plans 
to attempt the detection of these lines with ground-based near-IR observations (searching 
for the ~ 2 /im H2 lines) and possibly with the Herschel Space Telescope as a follow-up 
of our OT-1 imaging proposal. L Band VLA 21 cm observations could also be attempted 
in search of HI assuming that any structured emission around the star can be successfully 
disentangled from Galactic foreground and background emission. 



5. The Nebulosity Near the Two Stars 

Figure m shows the 8.0 and 24 /im emission in the region immediately around 6 Cephei 
and HD 213307. This is the area where the diffuse circumstellar emission is stronger at 
these two wavelength (box 1 in Figured]), and lowest at 70 /im, corresponding to a color 
temperature of ~ 180 K. The figure shows the rather complex morphology of this emission, 
with a small arched structure that seems to be connecting the two stars, and other filaments 
arching away from 6 Cephei. Can the morphology of this emission be explained in context 
of the 6 Cephei mass loss hypothesis? In absence of a complete geometry of the region (in 
particular, the relative position of the two stars along the line of sight) and kinematic data 
for the nebulosity, we can only speculate about the origin of these structures. 

One possibility is that the arched structure is in fact a second bow shock, this time 
associated to HD 213307. While the structure is not perfectly ahgned with the relative 
motion of HD 213307 with respect to the ISM (but still within the uncertainties), it is 
instructive to evaluate what would be the mass loss rate required for its creation. The 
standoff distance of the arc from HD 213307 is in this case dg ~ 26 arcsec (0.034 pc), and 
the escape velocity from HD 213307 (a B7-B8 main sequence star) is ~ 500 km s~^. Using 
equation [2l one obtains a required mass loss rate of 10"^° M© yr ^, more than one order 
of magnitude smaller than the minimum mass loss rate that we have estimated for 6 Cephei. 
This is in agreement wi th the expectation that a l ate-B star like HD 213307 should not have 



a strong wind (see e.g. iHempel &: Holweged 120031). as wind s drop off strongly in late B stars 
kudritzki fc PulslboOoT see also eq. 25 in lVink et al.ll200ll predicting log^^Q M ^ —11.6 for a 



B7 main sequence star with T^ff ~ 12, 500 K, L ~ 160 Lq, M ~ 4 Mq, solar metallicity and 
a galactic value of Voo/vesc = 1-3). The imperfect alignment of the arc with respect to the 
HD 213307 motion, and the brightening of the arc near 6 Cephei could be explained by the 
interaction of the two stellar winds and with the irradiation of the HD 213307 bow shock 
from the Cepheid photons. This explanation requires that the distance of the two stars along 



- 13 - 



the line of sight is significantly larger than their distance on the plane of the sky, outside 
the 3-dimensional paraboloid partially evacuated by the 6 Cephei wind. Conversely, if the 
two stars were on the same plane, a bow shock from HD 213307 would be blown away by a 
stronger wind from 6 Cephei. The uncertainly in the parallax of the two stars is large enough 
to allow this geometry. If this hypothesis is correct, then the possibility that 6 Cephei is the 
origin of a wind responsible for the large, 70 /xm bright, arc is strongly supported, as the 
weak wind from HD 213307 would not be sufficient to generate a bow shock at the standoff 
distance of 0.103 pc. 

Another possibility involves wind mass transfer between 6 Cephei and HD 213307, 
similar to the structu res observed between o Ceti and i ts compact corapanio n in the UV 



(iKarovska et al.l 119971 ) but also in the thermal infrared (iMarengo et al.ll200ll ). The much 
larger separation of 6 Cephei from HD 213307, with respect to the separation of o Ceti 
from its companion, makes however this interpretation very unlikely. An asymmetric mass 
ejection from 6 Cephei, in the chance direction of HD 213307 is also very unlikely because 
it would require non-radial oscillations with such strength that has never been observed in a 
Cepheid, or strong irradiation of the 6 Cephei atmosphere by the HD 213307 photons which 
is however negligible. The possibility that the arched structure is in fact originating from 
the FO V companion of HD 213307, being irradiated by the nearby B star, is also to be 
discounted. While the flux from HD 213307 contributes to as much as ~ 10% to the pho- 
tospheric flux of the FO V companion, a plume from this star (orbiting HD 213307 with a 



period of ~ 1 yr, see [Benedict et al.ll2002l ) would generate a tightly wound spiral. Assuming 
again a wind from HD 213307 as fast as its escape velocity, the spacing in the plume spiral 
would be ~ 100 AU, much smaller than the wide arc observed in the images. 

The exact nature of the structures seen in proximity of the two stars is difficult to 
assess with the available data. However, the presence of enhanced nebulosity nearest to the 
two stars is strongly supportive of the hypothesis that the overall structures that we have 
detected around 6 Cephei and its companion are local, rather than a chance alignement with 
ISM background cirrus. 



6. Discussion and conclusions 



Our S Cephei Spitzer data provides a new and unique view of the recent and current mass 
loss history of a Classical Cepheid star. This system is special, in that the presence of a B 
spectral type companion provides enough UV photons to excite the circumstellar molecules, 
making it possible to r nap them at infrared wavelengths. In the sample of 29 Cepheids 
observed with Spitzer (IMarengo et al.ll2010[ ). 8 have hot companions, but HD 213307 is 
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significantly hotter th an all of the others, with the exception of S Mus that also shows signs 
of extended emission (iBarmby et al.ll201l[ ). The question arises as to whether 5 Cephei may 
be unique in having mass loss as the result of the presence of HD 213307. The answer is 
most likely negative, given that the wide separation of the companion is probably too large 
to affect the mass loss properties of the Cepheid. None of the 5 Cephei properties as Cepheid 
(period, luminosity, or amplitude) suggest it is unusually likely to have a mass loss episode. 
We can then expect the mass-loss properties of 5 Cephei to be representative of Classical 
Cepheids in general, with the only difference that the presence of HD 213307 and the peculiar 
motion of 5 Cephei with respect to the ISM makes any stellar ejecta easier to detect. 

The presence of local (i.e. circumstellar) structures resembling a bow shock is strong 
evidence that one or both stars in the system are currently losing mass. The 70 /xm emission 
is better aligned with respect to the 5 Cephei position and relative proper motion vector 
(with respect to the ISM), although the low resolution and sensitivity at this wavelength is 
not good enough to rule out alignment with the companion HD 213307. While the wind 
responsible for this structure may b e arising from either star (or both), the circumstellar shell 
found by near-IR interferometers (IMerand et al.ll2006l ) around 5 Cephei strongly supports 
the possibility that the Cepheid star is the origin of the wind. The weak wind expected 
by a late B star like HD 213307 is also unlikely to be able to form a bow shock with the 
observed standoff distance of the 70 /im structure. An independent confirmation of this 
expectation could be obtained by showing that the arched filament near the two stars is 
indeed a secondary bow shock generated by a weak wind coming from HD 213307. 



Kervella et al.l (120091 ) found extensive (~ 100 to 1000 AU) "warm" (~ 100 K) emission 
in Spitzer and ground-based mid-IR images of RS Pup and £ Car, which they attributed to 
dusty mass loss from these stars (as opposed to cold dust emission, ~ 40 K, in the large 
massive nebula around RS Pup, attributed to ISM material compressed by the Cepheid 
wind). Our images suggest that part of the material responsible for the extended emission 
around 5 Cephei is of interstellar origin. Dust grains at ~ 80 K (the color temperature 
of the arc) are too cold to explain th e obse rved 8^.0/82^ ratio. On the other hand, PAH 
emission as described by iDraine fc Lil (120071 ) provides a reasonable fit of the data. Given 
the typical abundances of Cepheid atmospheres (where C/0 < 1), if dust were to form, 
its composition would largely be characterized by 0-rich silicates and not by carbon-based 
PAH. Furthermore, if indeed 5 Cephei {T^jf — 5800 K) is currently losing mass, the wind 
is likely to be atomic and dust free, at least in proximity of the stellar photosphere (e.g. 
Glassgold fc Hugginsl Il983l ). This is in ag reement with the ab sence of significant infrared 



excess we measured close to the star (see 



Marengo et al.ll2010l ). and the fact that the cir- 



cumstellar shell found by iMerand et al. (l2006h would have an equilibrium temperature too 
high for astrophysical dust survival. 
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We can conclude that the structures we detect with Spitzer are hkely the consequence 
of the interaction of a 5 Cephei wind with the local ISM. The low dust content of this wind 
near the star implies a different driving mech anism than th e dust -driven wind commonly 



associated to evolved stars. As suggested by iMerand et al.l (120061 ) . as an explanation for 
the shell structures observed with interferometry, it is conceivable that this "dustless" wind 
from S Cephei is perhaps triggered by the pulsation and shocks crossing the atmosphere. 
As it expands, this wind interacts with the ISM, leading to the formation of the bow-shock 
oriented with the proper motion of the Cepheid. The dilution of the ISM by the Cepheid 
wind is in agreement with the observed low abundance of PAH (less than one third of the 
abundance in the NCG 7023 nebula). Without detailed dynamical data it is however difficult 
to accurately determine which fraction of the dust responsible for the strong 24 and 70 fim 
emission is of interstellar origin, and how much is condensed in the outflow as it expands 
and cools. 

A star with the evolutionary mass of S Cephei will stay in the instability strip for a 
tot al time trr ~ 1-5 M yr (second and third crossing of the strip during the star's blue loop. 



see 



Bono et al.l l2000l ). If mass loss with a rate at the upper limit of our measured range 
(M ^ 6 • 10~^ Mq yr~^) is sustained through all the Cepheid phase, the total mass that 
could be ejected would be as high as ~ 0.1 Mq. This is ~ 2% of the total current pulsational 
mass of 6 Cephei. While not enough to completely account for the Cepheid mass discrepancy, 
this value may represent a significant portion of it, with possibly more mass lost at the onset 
of central He-burning (RG) phase, before the star crossed the Cepheid instability strip. A 
mass loss rate closer to our lower limit (M ^ 5 ■ 10~^ Mq yr~^), however, would only result 
in the ejection of a more modest ~ 0.01 Mq, or ^ 0.2% of the current pulsational mass 
of 6 Cephei. Firm constraints on the Cepheid mass loss rate will also have a significant 
impact on the plausibility of the predicted Mass-Luminosity relation of intermediate-mass 
stars during central helium burning phases. 

Far-IR images with higher angular resolution, larger field of view and better sensitivity 
are required to provide a better estimate of the current 6 Cephei mass loss rate, and the 
contribution of mass loss to solve the Cepheids mass discrepancy. Targeted observations 
with the Herschel space telescope, and high sensitivity radio maps in atomic and molecular 
H and CO, may be capable of characterizing the morphology and thermal structure of the 
large scale emission around this star, thus probing its mass loss history before and during 
the onset of the Cepheid phase. 

In summary, our Spitzer images have detected strong extended nebulosity around 6 Cephei 
and its companion HD 213307. This nebulosity is likely to be local, due to mass loss pro- 
cesses from one or both stars. A large scale far-IR arc, in particular, may be associated with 
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a bow shock generated by the wind of 5 Cephei as it moves through the ISM. This is one 
of the strongest direct evidence to date of Cepheid mass loss, and the best available direct 
measurement of the mass loss rate for a Classical Cepheid, which we estimate to be in the 
range oi ^ 5 ■ 10'^ to 6 • 10"^ Mq yr'^ 

This work is based on observations made with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory, California Institute of Technology under NASA 
contract 1407. P.B. acknowledges research support through a Discovery Grant from the 
Natural Sciences and Engineering Research Council of Canada. NRE acknowledges support 
from the Chandra X-ray Center, NASA contract NAS8-03060. This work was supported 
in part by the NSF REU and Department of Defence ASSURE programs under grant no. 
0754568. 

FaciHties: Spitzer (IRAC, MIPS). 
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Table 1. Observations Log and Magnitudes 



JD-2,400,000 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


[24] 


[70] 


6 Cephei 
53935.018^ 
53957.101'= 
54732.332'= 


2.174±0.044 
2.372±0.014 


2.183±0.037 
2.392±0.015 


2.166±0.032 
2.358±0.019 


2.150±0.039 
2.345±0.033 


2.120±0.001 


2.30±0.59 



HD 213307 

53935.018^1 ••• ••• ••• ••• 6.139±0.087 

53957.10r 6.391±0.010 6.399±0.010 6.413±0.006 6.406±0.005 

54732.332'= 6.237±0.137 6.383±0.117 6.402±0.304 6.363±0.376 



^ 20" radius aperture, with a 40"- 50" sky annulus and aperture correction of 1.13 at 24 /im; 16" 
radius aperture with 18"- 39" sky annulus and 2.04 aperture correction at 70 /im. 

'= PSF- fitting photometry as described in Marengo et ah ( 2O10l) . 

6.1" radius aperture, with 40"-50" sky annuhis and aperture correction of 1.048, 1.052, 1.053 and 
1.062 at 3.6, 4.5, 5.8 and 8.0 respectively. The photometry has been corrected for the field dependent 
distortion, and pixel -phase dependent uncertainties were added at 3.6 and 4.5 /im (as explained in 
Marengo et all |2010l ). 



7" radius aperture, with 40"-50" sky annulus and 1.61 aperture correction. 



Table 2. Extended Emission Surface Brightness' 



Box 


Ss.e 


S4.5 


S5.S 


'S's.o 


^24 


5*70 (DAT default)*^ 


570 (DAT custom)'^ 




[MJy sr-i] 


[MJy sr-i] 


[MJy sr-i] 


[MJy sr"i] 


[MJy sr-i] 


[MJy sr^i] 


[MJy sr-i] 


1 


< 0.33 


< 0.27 


0.82 ±0.12 


3.56 ±0.20 


5.13 ±0.27 


2.38 ±0.72 


1.69 ± 1.02 


2 


< 0.18 


< 0.18 


0.45 ±0.07 


1.60 ±0.08 


2.53 ±0.14 


4.58 ±0.61 


4.52 ±0.80 


3 


< 0.42 


< 0.33 


< 0.36 


0.79 ±0.08 


1.11±0.11 


4.78 ±0.78 


3.77 ±0.92 


4 


< 0.42 


< 0.33 


< 0.36 


0.63 ±0.07 


1.29 ±0.11 


5.23 ±0.62 


8.41 ±0.91 


NGC 7023 


53.2 ±7.5 


33.0 ±3.4 


300. ± 24. 


769. ± 60. 


481. ±40. 


2840 ± 180 





^ Values without uncertainty are 3cr upper limits. 

^ From the DAT reduced image, standard filtering (panel b Figure [2]). 

^ From the DAT reduced image, custom filtering (panel c Figure [2]). 
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Fig. 1.— IRAC (2008 September 23) and MIPS (2007 July 19) images of 6 Cephei and 
surrounding area. All images are shown in a linear surface brightness scale ranging from to 
8 MJy sr-^ Both S Cephei (A) and HD 213307 (B) have been PSF subtracted in all images 
except MIPS 70 fim: their location is identified either by the core PSF subtraction residuals, 
or by a cross. All background stars have been subtracted in the 24 fim image. The large 
arrow shows the 6 Cephei proper motion direction relative to its local ISM. Diffuse emission 
surface brightness was measured in boxes 1 to 4, with sky level estimated from boxes si and 
s2. The total 70 /im flux within the arched structure was derived in the dotted rectangular 
aperture (large box, excluding the small box centered on the S Cephei location) plotted in 
the last panel. The dark and clear cross-bars, centered on 6 Cephei in the IRAC images, are 
column pulldown and banding artifacts. 
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Fig. 2. — MIPS 70 /iin 6 Cephei image processed with the DAT piphne: (a) no filtering (re- 
taining strip pattern noise); (b) defauh filtering (as in post-BCD data); (c) custom filtering. 
All images are shown with the same linear color scale, from to 8 MJy sr^^. 
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Fig. 3. — Flux ratios of selected regions of the extended emission around S Cephei and on 
one bright filament in the reflection nebula NGC 7023. This figure makes use of the average 
values of the 70 fim surface brightness listed in Table [21 with their average errors. The 
Ss.o/S2A flux intensity ratio is a measure of the PAH fraction in the circumstellar material, 
while the 82^/810 reflects the dust temperature, driven by the stellar flux stochastically 
heating the grains. 
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Fig. 4. — Details of the IRAC 8.0 /iin (a) and MIPS 24 /im {b) images showing the area 
in proximity of 6 Cephei and HD 213307. Both stars have been PSF subtracted in the two 
images: 6 Cephei is the star north of HD 213307. All background stars have been removed 
in the 24 fim image to highlight the diffuse emission. The images are shown in linear color 
scale from to 8 MJy/sr surface brightness. The two arrows, having a length of 40 arcsec, 
point to the direction of the relative motion of the two stars with respect to the ISM. 



